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Abstract 


It is well known that highly polar and hydrogen bonding mixtures pose a serious 
challenge to equations of state. In the present work, it has been shown that excellent 
correlation and predictions of the behaviour for complex systems can be achieved 
when Wong- Sandler mixing rules are associated with a proper equation of state. In 
the pre.sent work, vapour-liquid equilibrium (VLB) calculations have been carried out 
for 15 binary and 5 ternary systems and the systems are^ chosen such that each 
system contains at least one polar substance. It is found that, Peng, Robinson, Stryjek 
and Vera (PRSV) equation of state along with Wong-Sandler mixing rules predicts 
VLB data which is in good agreement with experimental data. The Wong-Sandler 
mixing rules depend on two parann'ters, the excess molar Helmholtz free energy and 
second virial coefficient. The excess molar Helmholtz free energy is calculated by 
Nonrandom, Two-Liquid (NRTL) model and the cross second virial coefficient is 
calculated from Tsonopoulos correlation. 
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Chapter 1 


INTRODUCTION 


In chemical and petroleum industries, the separation of components from a 
mixture is one of the important tasks. In the manufacture of chemical products, 
several operations like extraction, adsorption, distillation, leaching and absorption 
are associated with phase-contacting. For the design engineer whose interest lies in 
sizing the equipment, predicting the operating costs and designing control schemes for 
new and existing separation processes, a knowledge of the thermodynamic behaviour 
of the systems becomes invaluable. 

Of all the separation processes, distillation is probably the most widely used in 
chemical industries, where the coexisting phases are vapour and liquid. Hence for the 
design of distillation columns, it is necessary to have reliable Vapour-Liquid Equilib- 
rium ( VLE) data, for a large number of binary and for a few multicomponent systems 
the experimental VLE data can be found in the literature. However tlv number of 
possible combinations of components of interest in chemical process design always far 
exceed the number of systems for which experimental VLE data are available. This 
necessitates the designer either to experimentally determine or theoretically estimate 
the required VLE data. 

One of the most important tasks of a design engineer is to develop a mathematical 
model to predict the behaviour of the systems. Among the several techniques of phase 
equilibrium calculations the analytic combined methods (Gamma-Phi approach) and 
the analytic direct methods (Equation of state approach) are widely used. There are 
two classes of thermodynamic models for phase equihbrium calculations: equations 
of .state and activity coefhcient (free energy) models. The regions of applicability of 
these models are schematically illustrated in Figure 1.1, where the axes are density 
and complexity of the mixture represented by the excess free energy in a high density 
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Figure 1.1, Schematic diagram of the approximate range of applica- 
bility of various thermodynamic models. 

d'he virial equation of state can be used at low densities to describe any mixture 
provided all the virial r.iefhcients ..le available. Activity coefficient models can also 
be used to describe the behaviour of liquid of any complexity, but only below its 
critical temperature. Cubic equations of state, with the usual Van der Waals one 
fluid mixing rules, can be used at all densities but only for relatively simple mixtures 
(i.e., hydrocarbons or hydrocarbons and inorganic gases) with low to moderate excess 
free energies. 

Equations of state applicable to both vapour and liquid phases have obvious ad- 
vantages over activity coefficient models, especially when one is interested in large 
ranges of temperature and pressure including critical and super critical regions. How- 
ever, phase equilibrium calculations by this method are rather lengthy except for 
simple equations of state. Another problem with the use of the equations of state 
technique is the insufficiency of information about the exact form of equations of 
state for mixtures and inaccuracy of the existing mixing rules for mixtures consisting 
of polar and associated molecules of large molecular size and shape differences. Re- 
cently, Wong and Sandler (1992) have developed a new theoretically correct equation 
of state mixing model which is not only conceptually simple, but more accurate for 
wider ranges of mixtures than other existing mixing models. 

In the present work, Wong and Sandler’s mixing rule is clubbed with Peng, Robin- 
son, Stryjek and Vera (PRSV) equation of state (or modified Peng-Robinson equation 
of state) to predict VLE data for few systems through f — <t> approach. In Wong- 
Sandlcr mixing rule the excess molar Helmholtz free energy and the activity coeffi- 
cients are calculated by Non Random Two Liquid (NRTL) model. The cross second 
virial coefficients Bij are calculated by Tsonopoulos correlation (1974). 
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Chapter 2 


THERMODYNAMICS OF 
PHASE EQUILIBRIUM 


2.1 The General VLE Problems 


In genieral, VLE problems are classified into four classes. They are, 

BIJBL P ; Calculate yk and P, given Xk and T 
DEW P ; Calculate and P, given yk and T 
BUBL T : Calculate yk and T, given Xk and P 
DEW T : Calculate Xk and 7’, given yk and P 

Thus, for a system containing N components, one specifies either P or P and 
either the liquid phase or the vapour pheise composition, fixing 1 + (A — 1) or A 
phase rule variables, which are the exact number required by the phase i ule, for the 
specification of two coexisting vapour-liquid phases. The total number of variables, 
to be determined, are A + 2, which includes vapour and liquid volumes. Hence, one 
needs A + 2 equations, to solve these A -f 2 variables. The required equations are as 
follows : 

(1) A equilibrium relationships for A components. 

(2) The equation of state for vapour phase. 

(3) The equation of state for liquid phase. 


Since all the (A -|-2) variables are involveld implicitly in these equations, one must 
(■m|)loy iterative procedure to determine these variables. 
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2.2 Thermodynamic Relationships and Equations 


'file internal energy of a closed, homogeneous system can be expressed as 

dU = TdS-PdV (2.1) 

Here U is the internal energy, S is the entropy, and V is the volume. For a system 
to be at equilibrium, the internal energy should reach a minimum value at constant 
S and V. That is 

dUsy = 0 (2.2) 

The following expression can be used in the calculation of energy changes, if 
the enthalpy //, is defined as: 

H = U + PV (2.3) 

From equations (2.1) and (2.3), one can get 

dH = TdS + VdP (2.4) 

To use temperature explicit and volume explicit expressions in the calculation 
of energy the Helmholtz energy, A is defined as 

A = U -TS (2.5) 

Differentiating and substituting equation (2.5) into equation (2.1) yields 

dA = -SdT - PdV (2.6) 

Equation (2.1) can also be rearranged so that the calculation of energy is made 
using expressions explicit in temperature and pressure. The arrangement is as follows 

G = H-TS (2.7) 

Here (? is the Gibbs free energy. Differentiation and substitution yields 

dG = -SdT + VdP (2.8) 

For equilibrium at constant temperature and pressure the Gibbs free energy must be 
at a minimum. That is 

dGx^p = 0 (2-9) 

2.3 Chemical Potential 

Equations (2.1), (2.4), (2.6), (2.8) are valid for systems containing only one com- 
ponent. For multicomponent systems, the internal energy is also a function of the 
number of moles of each component present. 

U = IJ{S, V'',ui,n2, . . . ,n,,^) (2.10) 
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If \v(' (k'firie chemical j)c)tcnial as 


■' ■ (£)„,, ■ (£)„.„ ■ (i)„, ■ (£)„,, »■“> 

Subslituting the definition of chemical potential equations (2.1), (2.4), (2.6), and (2.8) 
reduce to 


dU = TdS - PdV + Y. 

(2.12) 

dll = TdS + VdP + Y 

(2.13) 

dA = -SdT- PdV + Y 

(2.14) 

dG = -SdT + VdP + Y 

(2.15) 


t 


To luive a better understanding of the physical significance of the chemical poten- 
tial, one can inspect the terms of equation (2.12). If one considers the temperature as 
a measure of thermal driving force and pressure as a mechanical driving force, then 
the use of the term chemical potential makes more sense. When a closed, heteroge- 
neous system is at equilibrium, these potentials, or intensive properties, are equal for 
each component in each phase. If a system of m components has vapour and liquid 
phases in equilibrium, the following relationships hold : 


bn 

II 

(2.16) 

pV pL 

(2.17) 

p.\ = Mf 

(2.18) 

A = M2 



Mm = Mm 

2.4 Fugacity and Fugacity Coefficient 


The practical convenience dictates that the thermodynamic relationships for a sys- 
l-ern must be in terms of experimentally measurable quantities. In the case of phase 
equilibrium, the temperature and pressure can be directly measured, where as the 
chemical potential can not be. To overcome this difficulty, one must express the 
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chemical potential in terms of experimentally measurable quantities. By differen- 
tiating e(|uation (2.7) and combining the result with the first and second laws of 
thermodynamics the following equations can be derived. 

dfii = —SidT -t- VidP (2.19) 


( 2 . 20 ) 


Here s,- is the molar entropy and u,- is the molar volume. 
In the case of ideal gas, 

RT 


Substituting equation (2.21) in equation (2.20) and integrating at constant tempera- 
ture yields 


RTln- 


- /i? = RTln-^ (2.22) 

wliere the supercript “ signifies a reference state for the quantity. Although 
we do not have a direct means of calculating the chemical potential we do have an 
expression for the change in the chemical potential of an ideal gas as the pressure 
changes from some reference pressure P° to the pressure P. Equation (2.22) relates 
the chemical potential to the experimental quantities — pressure and temperature. 

Unfortunately, equation (2.22) is valid only for an ideal gas. To extend the rela- 
tionship to nonideal gases, G. N. Lewis defined a function /, called the fugacity, by 
the relation 

..0 nrnl^ /• ro r) 0 \ 




RTln- 


(2.23) 


FUr ideal gases, the fugacity is equal to the pressure, and for mixtures of ideal gases, 
the fugacity of a component is equal to its partiaU pressure. The fugacity is often 
referred to as a corrected pressure. If one considers two phases a and /? in equilibrium 
and sets 

/<“■ = f*"' (2.24) 

it follows that 


/r = f- 


(2.25) 


/“ = /f 


(2.26) 


The fugacity coefficient of a component i, (j)i is defined as 

V'P 


(2.27) 


In an ideal solution the fugacity coefficient of a component is equal to unity. The 
grcatc'r the fugacity coefficient deviates from unity, the greater the system deviates 
from ideal solution behaviour. 
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2.5 Equation of State Approach 


Starting with the definitions of enthalpy, entropy, Helmholtz energy and Gibbs free 
energy, we can arrive at the following relationship for the fugacity of component i in 
a mixture : 


RTlnc^i = RTln 


A 

ViP 


/ 

Jo L 


RT 

P 


dP 


Here, Vi is the partial molar volume of i which is defined 

- 

V / T,P,nj 


U, = 


(2.28) 


(2.29) 


For pure components, the partial molar volume u,- is equal to the molar volume u,, 
and equation (2.28) reduces to 


f 

= /„ 



(2.30) 


Unfortunately, explicit representation of volume by an equation of state is often dif- 
ficult. A better alternative is to characterize the system with an equation of state 
explicit in pressure. For a given volume and temperature, the system can exist at only 
one pressure. In contra.st, volume is not a unique function of temperature, pressure 
and concentration. For a pressure explicit equation of state, equation (2.28) can be 
written as 


RTln4>i = RTln^ = 
ViP 



dV - RTlnz 


(2.31) 


Here z is the compressibility factor defined as z = ^ 


Equation (2.26) requires that for vapour-liquid equilibrium the fugacity of a com- 
ponent in the vapour phase is equal to the fugacity of that component in liquid phase. 
Substituting the definition of the fugacity coefficient in equation (2.26), one gets 

(2.32) 


or = Xi4>f ’ (2.33) 

Here, y,- is the mole fraction of component i in the vapour phase and a:,- is the mole 
fraction of component i in the liquid phcise. 


In some cases the constants in the equation of state can be calculated solely from 
information about pure components constituting the system. Thus, these equations 
are predictive. Often the accuracy of the equation of state can be greatly enhanced 
by the introduction of some adjustable parameters, the values of which depend on 
the conditions of the system. With the introduction e'' adjustable parameters, the 
equation of state becomes a correlative method as opposed to a predictive one. In- 
formation about the system at discrete values of T, P and concentration is used to 
find optimum adjustable parameters, which in turn can be used for interpolation or 
extrapolation. 
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Chapter 3 


EQUATIONS OF STATE AND 
MIXING RULES 


By equations of state we refer to volumetric equations of state, that is, relationships 
between the pressure P, molar volume u, and absolute temperature T. From such a 
volumetric equation of state we get not only P — v — T information but, also (from 
the interrelationships provided by classicial thermodynamics) departure functions and 
phase equilibrium data. Despite the abundance of such equations of state, none has 
proven completely satisfactory for all areas of interest in scientific and engineering 
applications. 


3.1 Ideal Gas Equation 


The simplest volumetric equation of state for a gas, consistent with Boyle’s law and 
Charles’ law, is the ideal ga.s equation which is given by 


Pv = RT 

(3.1) 

Pv 

(3.2) 

or z = == 1 

RT 


where R is the universal gas constant and z is the compressibihty factor. While this 
equation of state is satisfactory for nonassociating gases at low pressures, it does not 
give an adequate description of the behaviour of gases at moderate and high pressures, 
and it does not predict the occurance of a vapour-liquid phase transition. 
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3.2 Cubic Equations of State 


!‘'or an accurate description of the P — v — T behaviour of fluids over wide ranges 
of temperature and pressure, an equation of state which is more comprehensive than 
tlie virial equation of state is required. Such an equation must be sufficiently general 
to apply to liquids as well as to ga.ses and vapours. It must not be so complex as 
to present excessive numerical or analytical difficulties in applications. Polynomial 
equations that are cubic in molar volume offer a simplicity, which are in fact the 
simplest equations capable of representing both liquid and vapour phases. 

But the simple equations are not accurate enough for the many P — v — T rep- 
resentations, and a more complex form of equation such as Ben edict- Webb- Rubin 
equation of state or Beattie- Bridgeman equation of state ( which are having eight 
and five adjustable parameters respectively ) is required to describe accurately the 
P — V — T behaviour of real fluids. However, simple equations of state have found 
extremely valuable applications in phase equilibrium predictions, as demonestrated 
by Soavc (1972), Peng and Robinson (1976) and Benmekki and Mansoori (1987) 
etc. 


3.2.1 Van der Waals Equation of State 


The first equation of state to give a qualitative description of the vapour and liquid 
phases and phase transitions was the famous cubic equation of Van der Waals (vdw), 
proposed in 1873. 

p = 7 -^ - 4 ( 3 - 3 ) 

[V — b) 


In this equation tlie constant b is the excluded volume, that is, that part of the 
molar volume which is not available to a molecule due to the presence of others. This 
contribution increases the pressure above that for an ideal ga^ at the same density 
and temperature. The second term on the right, which \fan der Waals took to be 
independent of temperature, is due to attractive forces and decreases the pressure. 

Van der Waals observed that the critical isotherm had a horizontal slope and an 
inflection point at the critical point, and utilized these characteristics in determining 
the parameters. Mathematically the conditions are expressed eis 



(3.4) 

(3.5) 
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(3.6) 


riuMi l)y applying these conditions to equation (3.3) one can get 



a = (27/64) R^TljP, 


(3.7) 


While Uan der Waals equation of state is of historical interest, it is not quanti- 
tatively accurate. For example, it predicts that the critical compressibility. 


Pc'^c 

RZ 


(3.8) 


is 0.375 for all fluids, while the value for different hydrocarbons varies from 0.24 to 
0.29; the range is wider when nonhydrocarbons are considered. Also, the predicted 
vapour pressures by tfcn der Waals equation of state are inaccurate. 


3.2.2 Redlich-Kwong (RK) Equation of State 


An important modification of the Van der Waals equation of state was made by 
Redlich and Kwong (1949), who introduced a temperature dependent term and a 
slightly different volume dependence in the attractive term. The RK equation of 
state is given by 


RT _ a/T°-^ 
(v — b) [u(u + 6)] 


( 3 , 9 ) 


It is belived to be the first simple equation of state that found successful industrial 
applications. The improvement wa.s so remarkable that many investigators have since 
attempted to make further improvement in the equation of state with varying degrees 
of success. This equation gives a somewhat better critical compressibility (z^ = 1/3), 
and better second virial coefficients {Tsonopoulos and Heidman, 1985) than the Van 
der Waals equation, but it is still not very accurate for the phase boundary (vapour 
pressure) and liquid density. 


3.2.3 Soave-Redlich-Kwong (SRK) Equation of State 


Wilson (1964) changed the temperature dependence of the attractive parameter in 
Redlich-Kwong equation by writing 


RT 


<T) 

b v{v -b b) 


( 3 . 10 ) 


where, a{T) = a(Tc).a(Tr,u;) 


(3.11) 
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a{T,) = 0.42747^ 

^ C 


(3.12) 


6 = 0.08664 


RT, 

Pr 


a = Tr 


l + (1.57 + 1.62u;)(;^-l 

1 r 


where u is the acentric factor {Pitzer, 1955, 1977) defined as 

u; = - 1.0 and Tr ^ 

However, it was Soave’s modification (1972) for the a parameter, 
Q = [1 + (0.48 + 1.57a; - 0.176a;^)(l - 


(3.13) 

(3.14) 


(3.15) 


(3.16) 


resulting in accurate vapour pressure predictions (especially above 1 bar) for light 
hydrocarbons, which led to cubic equations of state becoming an important tool for 
the prediction of vapour-liquid equilibria at moderate and high pressures for nonpolar 
fluids. This modification is referred to as the Soave-Redlich-Kwong or SRK equation. 


3.2.4 Peng-Robinson (PR) Equation of State 


There are some shortcomings which the Soave-Redlich-Kwong equation and the orig- 
inal Redlich-Kwong equation have in common. The most evident is the failure to 
generate satisfactory density values for the liquid even though the calculated vapour 
densities are generally acceptable. The other shortcoming is the prediction of an un- 
realistic universal critical compressibility factor of 1/3 for all substances. To alleviate 
these shortcomings, Peng and Robinson (1976) used a different volume dependence to 
give slightly improved liquid volumes (i.e., = 0.307) and changed the temperature 

dependence of a to give more accurate vapour pressure prediction for hydrocarbons 
in the 6- to 10- carbon number range (i.e., gcisoline fractions) as follows : 


p _ RT 0^0 

V — b v(v -f 6) -h b{v — b) 

with 


(3.17) 


Q = [1 + (0.37464 + 1.54226a; - 0.26992a;2)(l - r"-®)f (3.18) 


Uc 


0.45724 


R?T^ 

Pc 


(3.19) 


b = 0.07780^ 


(3.20) 
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3.2.5 Peng- Robinson-Stryjek- Vera (PRSV) Equation of State 


Although in many respects the modifications introduced in the modified Pcng-Rohinson 
(PRSV) equation of state by Stryjek and Vera (1986) follow ideas of previous workers 
in the field. The differences in the details are significant enough to produce a defi- 
nite improvement with respect to other versions of cubic equations of state. In this 
equation all other parameters are the same as the Peng-Robinson equation of state 
except for a. The parameter a is given by 

a=[l+k{l- r°-®)f (3.21) 

where k is considered to be a function of acentric factor cu only. After a careful 
examination of the deviations in the calculated vapour pressures at low reduced tem- 
peratures for compounds with a wide range of acentric factors, they have modified 
the functional dependence of k as 

k = ko + k,{l + T°-^){Q.7 -Tr) (3.22) 

with 

ko = 0.378893 -f 1.4897153cu - 0.17131848a;2 + 0.0196554a;^ (3.23) 

and ki, being an adjustable parameter characteristic of each pure substance. 

The Peng-Robinson (PR) and SRK equations are widely used in industry, espe- 
cially for refinery and reservoir simulation. The advantages of these equations are 
that they require little input information (only the critical properties and acentric 
factor for the generalized parameters), little computer time, and, for hydrocarbons, 
lead to good phase equilibrium predictions which is important for process design. 
However, these equations do have some important shortcomings. For example, liq- 
uid densities are not well predicted, the generalized parameters are not accurate for 
wonhydrocarbons (especially polar and associating fluids), and these equations do not 
lead to accurate predictions for long chain molecules. Also these equations are not 
accurate in the critical region, and vapour pressure predictions are not very accurate 
below 10 Torr. Even so, the Soave and Peng-Robinson equations are the most widely 
used cubic equations of state in industry. 


3.3 Virial Equation of State 


All the foregoing equations of state are empirical or at best semi- theoretical, but 
they are in closed form, which is of great advantage for practical applications. In 
contrast to the equations of state mentioned earlier, the virial equation of state has 
a theoretical basis from statistical mechanics, but has the form of an infinite power 
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HC'ries in volume, which is not convenient for practical use. 
state is given by 


P = RT 



B C 
^ H — ;• + 



The virial equation of 


(3.24) 


'Fhe virial coefficients B, C, D etc. are dependent only on temperature. This 
equation truncated after the second or third term depending on the availability of 
the virial coefficients is occasionally used for vapour pha.se properties calculations. 
Although, not widely used for predicting the phase equilibrium data, this equation is 
worth reviewing, because all the foregoing equations bear a formal resemblance to it 
when expanded into infinite series. For instance, if Van der Waals equation of state 
is expanded into infinite series, we get 


P 


RT 


-2 = 

yl 


1 ^ — 7 ^ }P 

- + ^ + -r + ... 


(3.25) 


All other closed forms of equations of state can be similarly expanded. The virial 
equation, originally suggested by Onnes as early as 1901, can be derived by the 
method of statistical mechanics. These methods show that the virial coefficients are 
precisely related to the molecular interactions of particular number of molecules. For 
instance, the second virial coefficient reflects the binary interactions and the third 
virial coefficient reflects the triple interactions and so on. 


3.4 Mixing Rules for Cubic Equations of State 


The equations of state are generally developed for pure fluids first, then extended to 
mixtures. The extension requires the so called mixing rules, which are simply means 
of calculating mixture parameters equivalent to those of pure substances. So far no 
exact theory like that for the virial equation prescribes this composition dependence, 
and quite often one relies on empirical mixing rules to provide approximate relation- 
ships. Recently Wong and Sandler (1992) have developed a new theoretically correct 
equation of state mixing model which is more accurate for a wider range of mixtures 
than other mixing models available. 


3.4.1 Van der Waals Mixing Rules 

The conventional Van der Waals type mixing rules are given by 

O-m — XjXjajj 

i i 

and bjn = 

» j 


(3.26) 

(3.27) 
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aij - y/a^j{l - kij) 


(3.28) 


% = ' (3.29) 

where kij = kji and hij = hji are binary interaction parameters to be determined 
from the experimental vapour-liquid equilibrium data. 


3.4.2 Linear Mixing Rules 


The livid r mixing rules are given by 

(Ljn ^ ^ ^ ^ X^Xj(l{j , dij yJCL^CLj 

t J 


(3.30) 


and 


i 


(3.31) 


3.4.3 Cubic Mixing Rules 


The cubic mixing rules are given by 


a 


1/3 ^ 
m 


a:i(aiiMi)^/^ + X2(a22M2)^''^ 


+ [(aiia22)^^^]^''^a;iX2(l - A:i2) 


(3.32) 


2 2 

and bm = (3.33) 

.'=1 j=i 

zvhere, bij = -Ct — _ fi-j'j (3.34) 

kij and hij are binary interaction parameters to be determined from experimental 
data. Ml, M 2 are molecular weights. 


3.4.4 Wong-Sandler Mixing Rules 

The Wong-Sandler mixing rules have been described in detail in chapter-4. 
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Chapter 4 

WONG- SANDLER MIXING 
RULES FOR PRSV EQUATION 
OF STATE 


While simple cubic equations of state can describe the behaviour of pure fluids 
reasonably well, they give satisfactory descriptions only for relatively simple mix- 
tures. It has long been recognized that the source of this difficulty is the mixing and 
combining rules which are used to obtain the equation of state iiuxture parameters 
from the pure fluid parameters. Consequently, much effort in recent years has been 
devoted towards developing alternate mixing rules. 

Several authors (Panagiotopoulos and Reid, 1986; Adachi and Sugie, 1986; San- 
doval et ai, 1989; Schwartzentruber et al., 1989) have proposed various modified 
forms of the Van der Waals mixing rules that use composition-dependent binary in- 
teraction parameters. For some highly non-ideal mixtures these rules have been used 
successfully. However, they are not generally applicable because in the low density 
limit they are inconsistent with the statistical mechanical result that the second virial 
coefficient must be a quadratic function of composition. Further, these rules may fail 
for simple mixtures {Shibaia and Sandler, 1988). 

To overcome these difficulties attempts have been made to develop density-dependent 
mixing rules (for example, Luedecke and Prausnitz, 1985; Panagiotopoulos and Reid, 
1986) so that the correct low-density limit is recovered. Such an approach, however, is 
adhoc and does not preserve the cubic nature of equation of state when used for mix- 
tures. There have been ongoing efforts to correct such deficiencies {Schwartzentruber 
and Renon, 1991; Mathias et ai, 1991), but with little success. 

.'\n alternate approach for developing mixing rules was first proposed by Huron 
and Vidal (1979). This method involves equating the excess Gibbs free energy at 
infinite pressure calculated from an equation of state to that calculated from the 
excess Gibbs free energy (activity coefficient) models for liquids. In order that the 
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('X('('ss (libbs ire(^ energy at infinite pressure computed from an equation of state be 
finit(q tluron and Vidal had equated the excess volume at infinite pressure to zero; 
otherwise, the Pv^ term becomes infinite. Consequently, a linear mixing rule for the 
excluded volume parameter 6^ had to be used. It was pointed out by Sheng et al. 
(1992) that such an assumption is not necessary if one equates excess Helmholtz free 
energy instead of excess Gibbs free energy . 

The Huron and Vidnl mixing rule does not satisfy the requirement that the second 
virial coefficient is a quadratic function of composition, and therefore, is inconsistent 
with statistical mechanical theory. Wong, Orbey and Sandier (1992) have shown from 
experimental data that 

g''‘ (7’, 7^=1 bar, x) = a^ (7’, P=\ bar, x) = a® (T, P=high pressure) 
riiis provides the justification for using exactly same functional form for at infinite 
pressure as has been used for g^ at low pressures. Further, it also explains why in 
Wong Sandler’s mixing rule model parameters obtained from low pressures data can 
be used, without change, at high pressures in a^. Consequently, this new class of 
mixing rules provide results that are correct at both the low density and high density 
limits without being density dependent. Moreover, these mixing rules show the correct 
(luadratic composition dependence of the second virial coefficient. 


4.1 The Mixing Model 


It is important to accurately estimate the pure component vapour pressures for pre- 
dicting the vapour-liquid equilibrium data with accuracy. Many of the components 
considered in this work do not fit the generalized correlation for the equation of 
state parameters of Peng and Robinson. Therefore, the present work considers the 
correlation of Stryjek and Vera (1986) which is given below. 


V — b v(v -b f>) + b(v — 6) 


(4.1) 


In this equation of state, P, T, v and R are the pressure, temperature, molar volume 
and universal gas constant respectively. The parameters a and b are given by 

a = (0.457235P"T,V Pc)a(T) (4.2) 


b = 0.077796PTe/Pc 


(4.3) 


a{T) = [l -b k{l - r°-®) 


(4.4) 


k = ko + k,{l+T^-^){0.7-Tr) 


(4.5) 


and, A:o = 0.378893 -bl.4897153w- 0.17131848u;2 + 0.0196554u;^ (4.6) 
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wlu'r(‘ uj is the acentric factor, is the critical temperature, Pc is the critical pressure, 
/ and ki is a constant specific to each fluid. 

If we exp<iiid P/LS equation of state as a virial series, the relation between the 
second virial coefficient B[T) and the equation of state parameters is given by 

= ( 4 . 7 ) 

hvom -Statistical mechanics^ the composition dependence of the second virial coefficient 
can be expressed as 


lWn^T.ldxiX,B.i(T) 

^ j 

(4.8) 


(4.9) 


wlu'ix' [() — a/ {ll'r)]ij is the composition independent cross second virial coefficient of 
th<' ecjuation ol state. A sufficient but not necessary, condition for satisfying equa- 
tion (4.9) is the {fav. dcr Waals one fluid model. Traditionally this solution has been 
used at all densities. 


Huron and Vidal (1979) developed a mixing rule by requiring that the excess 
(!ibl)s free energy {(/'^) at infinite pressure (or v = b) computed from an equation of 
stat(' b(' equal to that computed from an activity coefficient model resulting in 


a = b 


- <7g^ 


( 4 . 10 ) 


where cr is a numerical constant which depends on the particular equation of state 
used. While this model, and its variations, have been successfully used for some 
complex mixtures, it is no better than the Van der Waals one fluid mixing rule 
for hydrocarbon systems. Further, Huron and Vidal mixing rules suffer from the 
drawback that they do not satisfy the second virial coefficient boundary condition. 


Recently Wowj and Sandler (1992) have shown that it is possible to satisfy the 
second virial coefficient boundary condition. Moreover they have also shown that the 
equation of state predicts the same excess Helmholtz free energy at infinite pressure 
a^ as a function of composition as is obtained from an activity coefficient model with 
the following mixing model : 


and 


Gjn 


= b 


m 


E 




(4.12) 
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wIk'ic' fx, (l(.‘|)(‘nd(’s on the equation of state used. The value of c is equal to [ln[ 2 ^'^ — 
1))/ 2" ’ ior PliSV ('(luation ol statin d'he cross second virial coefficient is given by 



( 4 . 13 ) 


where kij is now a binary interaction parameter for the second virial coefficient. 
Wong, Orbey and Sandler (1992) have shown that for highly nonideal mixtures, once 
this second virial coefficient binary interaction parameters are obtained at a single 
temperature or pressure, the results can be extrapolated with very good accuracy 
over a range ol 200^ C in temperature and hundreds of bars in pressure. 


4.2 Second Virial Coefficient Binary Interaction 
Parameter 


I'/ahank, Skyu and Ilanif {1995) have studied in detail the P-v-T behaviour of different 
polar and non[)olar systems and they have observed that for nonpolar systems 
is not greatly different from whereas for polar systems is greatly different 

from 'I'liey then realized that the approximation 

( 4 . 14 ) 

is satisfact.ory even when differs significantly from B"^^. Further, from experi- 

mental study of various systems they have seen that B^f^ equals to Bl^^ near Tr^i 2 
of 0.84 ± 0.04, where 

Tr.u = 7^ { 4 . 15 ) 

rcn2 = - K 12 ) ( 4 . 16 ) 

One can choose the value of this Aha as prescribed by Tsonopoulos, 1974 and Tsonopou- 
I0.S, 1979 or one can tune this value in order to match with experimental VLE data. 


In the present work, jBiP has been calculated by Tsonopoulos correlation. Tsonopou- 
los correlation for 8 x 2 in the case of nonpolar molecules is given by 




and for polar molecules it is given by 


j 5 (o) ^ q. 5 ( 2 ) 

RTr.. 


(4.17) 


( 4 . 18 ) 


where, 


= 0.1445 


0.330 0.1385 0.0121 0.000607 

7^8 

Is. 


rp rp2 


T’S 

fj, • 


( 4 . 19 ) 


18 



(4.20) 


/](» 


= 0.0637 + 


0.331 


0.423 


0.008 


7’8 


and 


= 


7^6 ^ 

r>u -^r^ix 


Neither a nor 6 can be estimated with much accuracy. Tsonopoulos^ however has 
correlated a and b for several compound classes, b is zero for nonhydrogen bonded 
materials, and a is given by the values or expressions. These expressions ; • c' functions 
of reduced dipole moment, which is defined as 


where 


Pr 


rp2 

c 


fi = dipole moment, debyes 
= critical pressure, bars 
1\. = critical temperature, K 


( 4 . 21 ) 


.Mi.xiiig ruloH for critical properties as prescribed by Tsonopoulos (1974) are as given 
l)(*low. 


7;p = //’ciTc.d - Kij) 

(4.22) 

n TTcij{PciVci/Tci + PcjVcjlT^j) 

(4.23) 

uji LOj 

- 2 

(4.24) 


'File following steps can be adopted to calculate the second virial coefficient binary 
int.eraction parameter. 

(1) Sot equal to 0.84 ±0.04. 

(2) Calculate 7^,12, Pc.,n and u;i 2 by Tsonopoulos mixing rules. 

(3) Find cross temperature, cT = Tr,u *Tc,i 2 

(4) I’ind and from equation of state at this cross temperature. Calculate 

by Tsonopoulos correlation at this cross temperature. Since, for almost all 
systems this cross temperature is very high, the polar contribution in Tsonopoulos 
correlation (i.e., B^'^'>) at this high temperature is negligible. This means that at higher 
temperatures polar molecules behave like nonpolar molecules. So, for correlation of 
we can use successfully Tsonopoulos correlation of nonpolar molecules for both 
polar and nonpolar molecules. 

( 5 ) = Btr 

(6) Calculate second virial coefficient binary interaction parameters as 

= 1 - (4.25) 

This kn t.an be used by extrapolating with very good accuracy over 200® C in tern- 
perature and 100 bar in pressure. 
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4.3 Liquid and Vapour Volumes 


I'he FliSV equation of state can be rewritten as, 

z^ + {B- l)z^ -{-{A- - 2B)z + {B^ + B^ - AB) = 0 

where, 

aP 

J12T2 

B = ^ 

RT 

Pv 
RT 

Tlie first step in the application of these equations is the calculation of liquid and 
vapour densities which satisfy the above equation for a given temperature, pressure 
and sets of vapour and liquid compositions. 

d’liis cubic equation in 2 has either three real roots or one real root and two 
conjugate imaginary roots (Sokolnikoff, 1941). For the case where three real roots 
exist, the smallest root corresponds to the liquid phase and is called the liquid root 
and the largest root corresponds to the vapour phase and is called the vapour root.'^ 
The intermediate root has no physical meaning. As far as the mixtures are concerned, 
liquid root has been taken while solving for liquid phase and vapour root has been 
taken while solving for vapour phase. In case of one real root, it has been used as 
liquid root for Iquid phase and as vapour root for vapour phase. Newton - Rapkson 
iterative procedure has been used for solving the cubic equation. Iteration can be 
started from 2 = 0 for liquid phase and from 2 = 1 for vapour phase. 


(4.26) 

(4.27) 

(4.28) 

(4.29) 


4.4 Nonrandom, Two-Liquid (NRTL) Method 


The concept of local composition has been successfully used by Renon in his derivation 
of the nonrandom, two-liquid (NRTL) equation. However Renon’s equation, unlike 
Wilson’s, is applicable to partially miscible as well as completely miscible systems. 
For a solutioHof m components, the NRTL equation for the excess Gibbs free energy 
is given by 


EJL, r,,Gjai 
RT FA EteiGi.-ii 


(4.30) 


where, 


9ji 9ii 
RT 


(4.31) 
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(4.32) 


Gj. = exp {-OjiTji) (aji = aij) 


'I'lie activity coefficient for a component i is given by 

_ ET=1 

UUGux^ Er=iGo-^/ ) 


(4.33) 


(jij is an energy parameter characteristic of the i-j interaction. The parameter 
aij is related to the nonrandomness in the mixture. When aij is zero the mixture is 
completely random and equation (4.30) reduces to the two suffix Margules equation. 
For moderately nonideal systems, the NRTL equation offers no advantages over the 
simpler van Laar and three suffix Margules equations. However for strongly nonideal 
mixtures, and especially for partially immiscible systems, the NRTL equation often 
provides a good representation of experimental data provided care is exercised in data 
reduction to obtain the <Kljustable parameters. 


4.5 Other Necessary Expressions For VLE Cal- 
culation 


The parameter a, appearing in equation (4.10) is a numerical constant which depends 
on the particular equation of state used. The value of cr for PRSV equation of state 
is given by 

a = -^ln(\/2-l) (4.34) 

\/2 

The excels molar Helmholtz free energy at infinite pressure is given by 


aRT 


bmRT 


biRT 


(4.35) 


The above expression can be split into the following two equations 

(1 - D) 

with Q and D defined as : 

I J 

D - Yx - 

^ 'biRT aRT 

l * 


(4.36) 

(4.37) 


(4.38) 

(4.39) 
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When equation of state is used, equation (2.31) reduces to 


In <pi 


In 


P{v - b,n) 


RT 


+ 


j_ / dnbm \ f Pv _ W , ^ ( °-m 

b^[dn, \RT J 2V2\b^RT 


I fldn'^ 


1 f dnbr, 


In 


V + bmjl - V2) 

y + ^’m(l + \/2) 


\n drii j h^\ dui 
The partial derivatives of and bm are given by 

dnK„ 1 Q / dnD\ 

dn, ~ {I- d)[ti dui j (1 - oy V / 


and 


1 / 1 OrPa^, 


= D 


dnbr, 


RT \n dni ) dni 
The partial derivatives of Q and D are given by : 


+ br 


dnD 

^~d^ 


1 drPQ ' 
n dn; 




dnD 


dn, biRT a 


a,- _^ln 7 ?‘ 


with 


ln 7 r 


1 dna^ 
RT dn. 


(4.40) 

(4.41) 

(4.42) 

(4.43) 

(4.44) 

(4.45) 
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4.6 Block Diagrams 
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Figure 4. 1 Algorithm for bubble temperature calculation 
using an Equation Of State 



























Chapter 5 


RESULTS AND DISCUSSION 


In the application of the analytic direct method for vapour-liquid equilibrium 
calculation, the success of the vapour-liquid equilibrium predictions depend on the 
accuracy of the equation of state and on the mixing rules that are used. The PRSV 
equation of state predicts vapour pressure of a component i^aiXre accurate.' than 
other existing cubic equations of state. The Wong-Sandler mixing rule is theoretically 
correct which is more accurate for a wider range of mixtures than other mixing models 
available. In the present work, the Wong-Sandler mixing rule along with the PRSV 
equation of state is used to predict VLE data for few binary and ternary systems 
through (j) — <f> approach. In Wong-Sandler mixing rule, the excess molar Helmholtz 
free energy and the activity coefficients are calculated by NRTL model just because 
of its simplicity. 

In the following pages, the results of 15 binary systems and 5 ternary systems 
of highly polar substances are presented in tabular form as well as graphical form. 
To test the applicability of this model to multicomponent systems, calculations have 
been [)erformed for 5 ternary systems. The systems are so chosen such that each 
system contains at least one polar substance. The cross second virial coefficients Bij 
are calculated by Tsonopoulos correlation in which binary interaction parameters for 
cross temperatures have been tuned in order to match calculated results with exper- 
imental results. 

For a total of 20 systems, both BUBL P and BUBL T computations have been 
carried out. For BUBL P the computed equilibrium pressures have been compared 
with the experimental equilibrium pressures and for BUBL T the computed equilib- 
rium temperatures have been compared with the experimental equilibrium temper- 
atures and the deviations have been calculated. Similarly, the equilibrium vapour 
compositions have been compared with experimental equilibrium vapour composi- 
tions and the deviations in the vapour phase compositions have been determined. 
The experimental data have been taken from ’’Vapour liquid equilibrium data col- 
lection”, DECHMA series and the data for critical properties have been taken from 
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”'rii(‘ Pr<)|>('rt,i(>s ol' Cjw's H/. I/K|ui<ls”. 


Priiusnit./, 


5.1 Tables 


ACETONE — BENZENE 
PIIESSUIIE OF THE SYSTEM (bar) 
1.0130 

CHANGED TSONOPOULOS-K VALUES 
.0000 .0000 


.0000 .OQpO 

TEMPCC) X 


t (exp) 

t (cal) 

xl 

x2 

80.1000 

80.6000 

.0000 

1.0000 

79.2000 

79.7000 

.0100 

.9900 

76.3500 

76.8000 

.0500 

.9500 

73.6000 

74.1000 

.1000 

.9000 

69.7000 

70.2000 

.2000 

.8000 

66.7500 

67.2000 

.3000 

.7000 

64.5000 

65.0000 

.4000 

.6000 

62.6500 

63.1000 

.5000 

.5000 

61.0000 

61.4000 

.6000 

.4000 

59.6000 

60.0000 

.7000 

.3000 

58.3500 

58,8500 

.8000 

.2000 

57.2500 

57.7000 

.9000 

.1000 

56.7000 

57.2000 

.9500 

.0500 

56.2700 

56.7000 

.9900 

.0100 

56.1800 

56.6000 

1.0000 

.0000 


AVERAGE DEVIATION IN TEMPERATURE CC) 


. 463 

AVERAGE DEVIATION IN Y-VALUES 
DYl DY2 

.00404 .00404 


Y. 

, exp 

yi 

y2 

.0000 

1.0000 

.0352 

.9648 

. 1496 

.8504 

.2531 

.7469 

.4030 

.5970 

.5140 

.4860 

.6030 

.3970 

.6785 

.3215 

.7464 

.2536 

.8100 

.1900 

.8737 

.1263 

.9371 

.0629 

.9687 

.. .0313 

.9937 

.0063 

1.0000 

.0000 


Y . cal 


yi 

y2 

.0000 

1.0000 

.0305 

.9695 

. 1366 

.8634 

.2415 

.7585 

.3951 

.6049 

.5073 

.4927 

.5975 

.4025 

.6755 

.3245 

.7463 

.2537 

.8127 

.1873 

.8765 

.1235 

.9387 

.0613 

.9694 

.0306 

.9939 

.0061 

1.0000 

.0000 


Table 5.1.1: Comparison of experimental and calculated values of y,- and t from Wong- 
Sandler mixing rules for Acetone-Benzene at 1.013 bar pressure. 


20 



ACETONE --- CHLOROFORM 
PRESSURE OF THE SYSTEM (bar) 


1.0130 

CHANGED TSONQPOULOS-K VALUES 


. .0000 

.0000 





.0000 

.0000 





TEMPCC) 

X 


Y. 

exp 

t ( exp ) 

t (cal) 

xl 

x2 

yi 

y2 

62.8000 

62.2000 

. 1108 

.8892 

,0650 

.9350 

63.1000 

62.5000 

.1375 

.8625 

,1000 

.9000 

63.8000 

63.2000 

.2108 

.7892 

. 1760 

.8240 

64.4000 

63.8000 

.2660 

.7340 

.2370 

.7630 

63.9000 

63.3000 

.4771 

.5229 

.5170 

.4830 

62.8000 

62.2000 

.5750 

.4250 

.6480 

.3520 

61.6000 

61.0000 

.6633 

.3367 

.7505 

.2495 

60.4000 

59.8000 

.7388 

.2612 

.8235 

.1765 

59.4000 

58.8000 

.7955 

.2045 

.8688 

.1312 

58.8000 

58.3500 

.8590 

.1410 

.9165 

.0835 

57.5000 

57.2000 

.9145 

.0855 

.9522 

.-0478 

AVERAGE 

DEVIATION 

IN TEMPERATURE CO 




• 559 

AVERAGE DEVIATION IN Y-VALUES" 
DYl DY2 

.01G20 .01620 


Y . cal 


yi 

y2 

0888 

.9112 

1152 

.8848 

1954 

.8046 

2614 

.7386 

5241 

.4759 

6377 

.3623 

7310 

.2690 

8031 

.1969 

,8525 

.1475 

,9031 

.0969 

,9438 

.0562 


T). 1 . 2 : ( !()m|)ii.riw>ii ol ('Xix'i in 


and .-aKadaU'd vajiu-s of and i from tf'r™;;- 


Sandkr mixiiis i i.l« lor ArdimcAMomJo,-,,, at 1.013 bar i>rcss,irc. 




ACETONE — HEXANE 
PRESSURE OF THE SYSTEM (bar) 
1.0130 

CHANGED TSONOPOULOS-K VALUES 
.0000 .1550 

.1550 .0000 


TEMP(JC) X 


l;(nx{)) 

1; (cnl) 

xl 

x2 

68.7500 

68.7500 

.0000 

1.0000 

58.9100 

59.0000 

.0942 

.9058 

53.5100 

53.6000 

.1986 

.8014 

51.5700 

51.6000 

.2989 

.7011 

50.5400 

50.5700 

.3983 

.6017 

50.0000 

50.1000 

.4995 

.5005 

49.7800 

49.8000 

.5999 

.4001 

49.8400 

49.8400 

.6999 

.3001 

50.3100 

50.3500 

.8007 

.1993 

51.7600 

51.8000 

.9015 

.0985 

56.2500 

56.3000 

1.0000 

.0000 

AVERAGE 

DEVIATION 

IN TEMPERATURE 


.045 

AVERAGE DEVIATION IN Y-VALUES 
DYl '■ DY2 
.00808 .00808' 


Y. 

exp 

Y. 

cal 

yi 

y2 

yi 

y2 

.0000 

1.0000 

.0000 

1.0000 

.3350 

.6650 

.3337 

.6663 

.4830 

.5170 

.4853 

.5147 

.5400 

.4600 

.5562 

.4438 

.5840 

.4160 

.5958 

.4042 

.6100 

.3900 

.6208 

.3792 

.6160 

.3840 

.6395 

.3605 

.6540 

.3460 

.6607 

.3393 

.7000 

.3000 

.6978 

.3022 

.7950 

.2050 

.7809 

.2191 

1.0000 

.0000 

1.0000 

.0000 


TmJ,Ic. 5. 1 of (•xi,(Tinu‘..l,al a..ul ca,lci.la,tc(l values of ?/, and t IVom Wony 
Sandler mixing rules for AceUmc-Ilexane. at 1.013 bar pressure. 
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ACETONE — METHANOL 
PRESSURE OF THE SYSTEM (bar) 
1.0130 


CHANGED 

TSONQPOULOS 

-K VALUES 


.0000 

.0500 



.0500 

.0000 


•' 

tempCc) 

X 


t (exp) 

t (cal) 

xl 

x2 

62.8000 

62.9000 

.0600 

.9400 

60.8000 

60.8000 

. 1500 

.8500 

58.6000 

58.7000 

.2950 

.7050 

57.2000 

57.3000 

.3980 

. 6020 

56.3000 

56.3000 

.5250 

.4750 

55.8000 

55.8500 

.7320 

.2680 

55.7000 

■' 55.8000 

.8000 

.2000 

55.9000 

56.0000 

.8920 

.1080 


AVERAGE DEVIATION IN TEMPERATURE 
.069 

AVERAGE DEVIATION IN Y-VALUES 
DYl DY2 

.00563 .00563 


Y. 

exp 

Y. 

cal 

yi 

y2 

yi 

y2 

0970 

.9030 

.1127 

.8873 

2500 

.7500 

.2503 

.7497 

4200 

.5800 

.4181 

.5819 

5170 

.4830 

.5127 

.4873 

6250 

.3750 

.6115 

.3885 

7540 

.2460 

.7548 

.2452 

8000 

.2000 

.8040 

.1960 

8750 

.1250 

.8795 

.1205 


'I'abler,.!.^: Comparison of oxpcnim(M.I,al anti calculaUrd v^vlncs-til^raud t from Wony- 
Sandlcr inixing rules for Acetone- Methanol at 1.013 bar pressure. 
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A(J1';'1'UN1': TK'rilACHLUllUMETHANE 

PRESSURE OF THE SYSTEM (bar) 

1.0130 

CHANGED TSONOPGULOS-K VALUES 
. 0000 - .0000 
.0000 ■ .0000 


TEMPCC) '■ 

X 


Y. 

exp 

Y. 

cal 

t (exp) 

t (cal) 

xl 

x2 

yi 

y2 

yi 

y2 

70.8000 

70.9000 

.0590 

.9410 

,2025 

.7975 

.1753 

.8247 

68.7400 

68.8000 

.0870 

.9130 

.2710 

.7290 

.2369- 

.7631 

64.4500 

64.5000 

.1790 

.8210 

,4075 

.5925 

.3827 

.6173 

61.9100 

62.0000 

.2640 

.7360 

.4895 

.5105 

.4739 

.5261 

59.8300 

59.9000 

.3740 

.6260 

.5655 

.4345 

.5626 

.4374 

58.7400 

58.8000 

.4510 

.5490 

.6125 

.3875 

.6149 

.3851 

57.9400 

58.0000 

.5255 

.4745 

.6550 

--■.3450-— 

-^6619 

.3381 

57.1800 

57.3000 

.6165 

.3835 

.7065 

.2935 

.7179 

.2821 

56.6700 

56.8000 

.6960 

.3040 

.7560 

.2440 

.7678 

.2322 

56.3600 

56.4500 

.7620 

.2380 

.7985 

.2015 

.8112 

. 1888 

56 . 1500 

56.3000 

.8295 

.1705 

.8460 

.1540 

.8586 

. 1414 

56.0100 

56.2500 

.8950 

. 1050 

.8980 

. 1020 

.9086 

.0914 

56.0200 

56.0320 

.9140 

.0860 

.9150 

.0850 

.9240 

. 0760 

55.9900 

56.1500 

.9530 

.0470 

.9520 

.0480 

.9571 

. 0429 

AVERAGE 

deviation 

IN TEMPERATURE CC) 






.099 

AVERAGE DEVIATION IN Y-VALUES 
DYl DY2 

.01336 .01336 


T'ahlci 5. 1 .5: Ooitiparisoii orc^xperiinciilal and calculalcd values oi y, 

Sandler mixing rules for Aceione-Tetrachloromeihane at 1.013 bar pressure. 
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, CARBON DISULPHIDE — ACETONE 
PRESSURE OF THE SYSTEM (bar) 

1.0130 

CHANGED TSONOPOULOS-K VALUES 
.0000 .1000 
.1000 .0000 

TEMP CO X Y.exp 


t (exp) 

t(cal) 

xl 

x2 

yi 

y2 

54.0000 

54.0000 

.0190 

.9810 

.0832 

.9168 

51.4000 

51.4500 

.0476 

.9524 

.1850 

.8150 

46.6000 

46.6000 

. 1340 

.8660 

.3510 

.6490 

44.0000 

44.1000 

. 1958 

.8042 

.4430 

.5570 

■V 

41.4000 

41.5000 

.2912 

.7088 

.5275 

.4725 

40.3000 

40.4000 

.3798 

.6202 

.5740 

.4260 

39.8000 

39.9500 

. 4477 

.•5523 

.5980 

.4020 

39.3000 

39.4500 

.5360 

.4640 

.6270 

.3730 

39.1000 

39.2500 

.6530 

.3470 

.6610 

.3390 

39.3000 

39.4000 

.7894 

.2106 

.7050 

.2950 

39.6000 

39.7000 

.8023 

.1977 

.7230 

.2770 

40.5000 

40.5500 

.8795 

.1205 

.7600 

.2400 

43.5000 

43.5500 

.9683 

.0.317 

A 

.8860 

.1140 


AVERAGE DEVIATION IN TEMPERATURE Cc) 
.085, 

AVERAGE DEVIATION IN Y-VALUES 
DYl '■ DY2 

.00557 .00557 


Y.cal 


yi 

y2 

.0752 

.9248 

.1682 

.8318 

.3532 

.6468 

.4352 

.5648 

.5193 

.4807 

.5704 

.4296 

.5994 

.4006 

.6291 

.3709 

.6614 

.3386 

.7043 

.2957 

.7098 

.2902 

.7573 

.2427 

.8915 

.1085 


T.,hl,. 5 I «■ ( of cxpcrimonlal .and calcnl.atcd values of and t from W.ng- 



ETHYLENE OXIDE — ACETALDEHYDE 





PRESSURE 

; OF THE SYSTEM (bar) 






1.0130 







CHANGED 

TSOMOPGULQS-K VALUES 






.0000 

.0200 







.0200 

.0000 







tempCc) 

X 


Y. 

exp 

Y . cal 

t (exp) 

t (cal) 

xl 

x2 

yi 

y2 

yi 

y^ 

17.3000 

17.3000 

.2950 

.7050 

3900 

.6100 

3896 

6104 

16.0000 

16.0000 

.4100 

.5900 

5130 

.4870 

5135 

4865 

14.6000 

14.6000 

.5280 

.4720 

6360 

.3640 

6282 

3718 

13.5000 

13.5000 

.6400 

.3600 

7200 

.2800 

7272 

2728 

11.6000 

11.6000 

.8830 

.1170 

9180 

.0820 

9172 

,0828 

10.7000 

10.7000 

.9666 

.0334 

9792 

.0208 

9768 

,0232 

10.7000 

10.7000 

.9675 

.0325 

9778 

.0222 

,9775 

.0225 

10.6000 

■' 10.6000 

.9794 

.0206 

,9861 

.0139 

, 9857 

.0143 

10.6000 

10.6000 

.9852 

.0148 

,9908 

.0092 

.9898 

.0102 

10.6000 

10.6000 

. 9877 

.0123 

.9912 

. 0088 

.9915 

. 0085 

10.6000 

10.6000 

.9895 

.0105 

.9930 

.0070 

.9927 

.0073 

10.6000 

10.6000 

.9912 

.0088 

.9947 

.0053 

.9939 

.0061 

10 .'1000 

10.4000 

.9917 

.0083 

.9956 

.0044 

.9943 

.0057 

10.4000 

10.4000 

.9965 

.0035 

.9982 

.0018 

.9976 

.0024 

AVERAGE 

DEVIATION 

IN TEMPERATURE (C) 





.000 







AVERAGE DEVIATION 

IN Y-VALUES 




— 



DYl 

DY2 






.00172 

.00172 






.'I 1 . 

7- CoinDarisoii of ex|)erim<'iil.a.l and 

caleul 

altul va,hK's ol 

' yi and t 1 

rum Wong- 

Sandler inixing rules lor hthylenc 

oxide- Acetaldehyde at 1.013 b 

? 

ar pressure. 



ACETALDEHYDE — ACETIC ACID 
PRESSURE OF THE SYSTEM (bar) 


1.0240 

CHANGED TSONOPOULOS-K VALUES 


.0000 

.0000 







.0000 

.0000 







TEMP CC) 


X 

Y. 

exp 

Y. 

cal 

t (exp) 

t(cal) 


xl 

x2 yl 

y2 

yl 

y2 

84.0000 

82.0000 


,0900 

.9100 .6380 

.3620 

.7344 

. 2656 

68.0000 

63.5000 


,1830 

.8170 .8080 

.1920 

.8860 

.1140 

58.0000 

55.5000 


.2460 

.7540 .8790 

.1210 

.9250 

.0750 

50.0000 

49.0000 


,3090 

.6910 .9270 

.0730 

.9474 

.0526 

42.0000 

42.0000 


.4190 • 

.5810 .9660 

.0340 

.9681 

.0319 

AVERAGE 

DEVIATION 

IN 

TEMPERATURE CC) 




2.000 







AVERAGE 

DEVIATION 

IN 

Y-VALUES 




DYl 

DY2 







.04857 .04857 


Tabl<.r, 1 8- (V«npari»ouc.roxpori.nontala,idcakulalcdvalucsof!,,a..daron. Wong- 

r,a AerUUlrMAM. „ad at .,024 bar pressure. 



ACETALDiaiYDE --- 'BENZENE 
PRESSURE OF THE SYSTEM (bar) 
1.0130 

CHANGED TSONOPOULOS-K VALUES 
.0000 .0300 

.0300 .0000 


TEMPCC) X Y.exp Y.cal 


t (exp) 

t (cal) 

xl 

x2 

yi 


yi 

y2 

80.1000 

80.6500 

.0000 

1.0000 

.0000 

1.0000 

.0000 

1.0000 

71.0000 

71.0000 

.0500 

.9500 

.3130 

.6870 

.3231 

.6769 

63.0000 

62.7500 

. 1000 

.9000 

.5220 

.4780 

.4926 

.5074 

54.2000 

53.0000 

.2000 

.8000 

.6940 

.3060 

.6692 

.3308 

49.4000 

47.0000 

.3000 

.7000 

.7670 

.2330 

.7671 

.2329 

46.2000 

42.0000 

.4000 

.6000 

.8160 

. 1840 

.8341 

. 1659 

44.0000 

37.0000 

.5000 

.5000 

.8540 

. 1460 

.8842 

.1158 

42.0000 

33.0000 

.6000 

.4000 

.8870 

.1130 

.9219 

.0781 

39.8000 

30.0000 

.7000 

.3000 

.9160 

.0840 

.9502 

.0498 

37.7000 

26.3000 

.8000 

.2000 

.9440 

.0560 

.9719 

.0281 

36.0000 

23.5000 

.9000 

.1000 

.9680 

.0320 

.9880 

.0120 

20.8000 

21.5000 

1.0000 

.0000 1 

.0000 

.0000 

1.0000 

.0000 

AVERAGE 

DEVIATION 

IN TEMPERATURE Cb) 






4.916 

AVERAGE DEVIATION IN Y-VALUES 
DYl .. DY2 
.01915 .01915 


T'abk'. 5. 1 .9: (Joiiiparisoii of cxpcrimciiial and calculated values ol y,- and t from Wong- 
Sandlv.r inixiiig rules for AccMildehyde- Benzene at 1.013 bar pressure. 



IJIETIIYL ETllEll --- BENZENE 
PRESSURE OF THE SYSTEM (bar) 
1.0160 

CHANGED TSONOPOULOS-K VALUES 
.0000 .0000 
.0000 .0000 


tempCc) X 


t (exp) 

t (cal) 

xl 

•' x2 

80.3000 

80.3200 

.0000 1 

.0000 

67.5000 

67.5000 

.1755 

.8245 

60.6500 

60.6700 

.2825 

.7175 

55.9500 

55.9500 

.3691 

.6309 

52.0000 

52.0500 

.4520 

.5480 

50.0500 

50.0300 

.4991 

. 5.009 

48.5000 

48.4700 

.5391 

.4609 

47.8000 

. 47.8400 

.5490 

.4510 

46.1000 

46.1200 

.5987 

.4013 

44.3500 

44 . 3800 

.6520 

.3480 

42.3000 

42.2500 

.7177 

.2823 

38.1500 

38.2000 

.8623 

.1377 

34.7500 

34.8000 

1.0000 . 

.0000 

AVERAGE 

DEVIATION 

IN TEMPERATURE 

.029 



AVERAGE 

DEVIATION 

IN Y-VALUI' 

:s 

DYl 

DY2 



.00402 .00402 


Y. 

exp 

Y. 

cal 

yi 

y2 

yi 

y2 

.0000 

1.0000 

.0000 

1.0000 

.4601 

.5399 

.4712 

.5288 

.6295 

.3705 

.6328 

.3672 

.7269 

.2731 

.7257 

.2743 

.7967 

.2033 

.7930 

.2070 

.8288 

.1712 

.8243 

.1757 

.8525 

.1475 

.8476 

.1524 

.8580 

.1420 

.8533 

.1467 

.8831 

.1169 

.8780 

.1220 

.9063 

.0937 

.9011 

.0989 

.9307 

.0693 

.9257 

.0743 

.9718 

.0282 

.9682 

.0318 

1.0000 

.0000 

1.0000 

.0000 


(“O 


Table 5.1.10: Coinparisou 
Wong-Sajidler mixing rules 


of cxperimeiilal and calculated values of y, and t 
for Diethyl cther-Benzene at 1.016 bar pressure. 


from 



DIETHYL ETHER — - CHLOROFORM 
TEMPERATURE OF THE SYSTEM (K) 
290.0000 

CHANGED TSONOPOULOS-K VALUES 


.0000 

.0000 



.0000 

.0000 



PRES 

S . (mm) 

X 


P (exp).i 

P(cal) 

xl 

x2 

143.000 

146.517 

. 1000 

.9000 

143.000 

15S.016 

.2000 

.8000 

151.000 

170.068 

.3000 

.7000 

166.000 

191.833 

.4000 

.6000 

193.000 

219.648 

.5000- 

.5000 

230.000 

252.208 

.6000 

.4000 

273.000 

287.813 

.7000 

.3000 

317.000 

324.650 

.8000 

.2000 

360.000 

361.019 

.9000 

.1000 

397.000 

395.544 

1.0000 

.0000 


AVERAGE DEVIATION IN PRESSURE (mm Hg) 
13.42331 


AVERAGE DEVIATION IN Y-VALUES 
DYl UY2 

0.02583 0.02553 


Y. 

GXp 

Y. 

cal 

yi 

y2 

yi 

y2 

.0980 

.9020 

.1281 

.8718 

.2240 

.7760 

.2899 

.7100 

.3970 

.6030 

.4602 

.5395 

.5720 

.4280 

,.6146 

.3847 

,7410 

.2590 

.7395 

.2595 

.8520 

. 1480 

.8331 

- . 16S7-— 

.9190 

.0810 

.8999 

.0990 

.9590 

.0410 

.9465 

.0528 

.9830 

.0170 

.9785 

.0213 

1.0000 

Toooo"’ 

“T.oooo 

.0000 


'l'ii,l)l<' 5,1.11: (l()in|)ariH()ii of ('.\i){’riiii(Mi(.a.l aiul ('a,lc.ula.t.('cl vaiiu's of t/, atul p li'oiu 
mixing rulon lor Diethyl ether- ( 'hi oi'ojovni a.t. 29U.0 K tciupcratauc. 



DIETHYL ETHER — - CARBON DISULFIDE 


TEMPERATURE OF THE 
293.0000 

CHANGED TSQNOPOULOS 
.0000 .0000 

.0000 .0000 

SYSTEM (K) 

-K VALUES 






PRESS 

. (mm) 

X 


Y. 

GXp 

Y. 

cal 

P(exp) 

P(cai) 

xl 

x2 

yi 

y2 

yi 

y2 

366.000 

335.954 

.1500 

,8500 

,2760 

.7240 

.2738 

.7292 

393.000 

359.761 

.2500 

.7500 

. 3820 

._.....618J_ 

.3878 

.6158 

421.000 

386.673 

.4000 

.6000 

.5060 

.4940 

.5219 

.4818 

432.000 

400.916 

.5000 

.5000 

.5790 

.4210 

.6008 

.4025 

443.000 

418.475 

.6500 

.3500 

.6750 

. 3250 

.7152 

.2873 

445.000 

432.309 

.8000 

.2000 

.8020 

. 1980 . 

.8321 

.1695 

446.000 

439.636 

.9000 

.1000 

.8990 

. 1010 

.9140 

.0870 

441.000 445.464 1.0000 .0000 1 

AVERAGE DEVIATION IN PRESSURE (mm Hg) 

.0000 

.0000 

1.0000 

.0000 


21.88020 

AVERAGE DEVIATION IN Y-VALUES ' 
DYl DY2 

0.018822 0.014288 


'I'ablo 5.1.12: (loiiipari.soii of cx[)crimcnlal aiul calculated val 
Won.(i-San.(ll('r mixing rules lor Diethyl cthcr-Carbon disulfide 
tur('. 


lies of yi and p Iroin 
at 293.0 K tempera- 
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DIETHYL ETHER ETHYL ACETATE 

TEMPERATURE OF THE SYSTEM (K) 

303.0000 


CHANGED TSOMOPOULOS-K VALUES 


.0000 

0000 



.0000 

0000 



PRESS 

! . (mm) 

X 


P(GXp) 

P(cal) 

xl 

x2 

227.000 

233.558 

.2000 

, 8000 

281.000 

288.416 

.3000 

.7000 

333.000 

342.573 

.4000 

.6000 

388.500 

395.731 

.5000 

.5000 

441.000 

447.736 

.6000 

.4000 

494.000 

498.624 

.7000 

.3000 

545.000 

548.688 

.8000 

.2000 

597.000 

598.563 

.9000 

.1000 

650.000 

644.499 

1.0000 

.0000 


AVERAGE DEVIATION IN PRESSURE(mm Hg) 


5.87954 

AVERAGE DEVIATION IN Y-VALUES 
DYl DY2 

.00445 .00426 


Y. 

exp 

Y. 

cal 

yi 

y2 

yi 

y2 

.5758 

.4242 

.5749 

.4265 

.7012 

.2988 

.6974 

.3038 

.7859 

.2141 

.7803 

.2208 

.8465 

.1535 

.8400 

.1609 

.8918 

.1082 

.8853 

.1155 

.9268 

.0732 

.9212 

.0795 

.9550 

.0450 

.9509 

.0496 

.9787 

.0213 

.9766 

.0237 

1.0000 

.0000 

1.0075 

.0000 
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ETHYL CHLORIDE DIETHYL 

TEMPERATURE OF THE SYSTEM (K) 


298.0000 

CHANGED TSONOPOULOS-K VALUES 


.0000 

.0400 



.0400 

.0000 



PRESS 

. (mm) 

X 


P (nxp) 

P(cal) 

xl 

■ x2 

967.000 

963.354 

.4530 

.5470 

1004.000.' 

1006.177 

.5210 

.4790 

1075.000 

1098.921 

.7000 

.3000 

1150.000 

1179.886 

.9200 

.0800 

1162.000 

1200.012 

1.0000 

.0000 


AVERAGE DEVIATION IN PRESSURE (mm Hg) 
16.97869 

AVERAGE DEVIATION IN Y-VALUES 
DYl DY2 

.01909 .01945 


ETHER 


Y. 

exp 

Y. 

cal 

yi 

y2 

yi 

y2 

.6667 

. 3333 

.6669 

.3381 

.7088 

.2912 

, .7116 

.2932 

.8068 

. 1932 

.8124 

.1913 

.8346 

. 1654 

. 9389 

.0623 

1.0000 

.0000 

.9985 

.0000 


Table 5.1.14: (loinparisoii 
Wong-Sandkr mixing rules 


ol expcriinculal and calculated values of yi and p from 
iov Ethyl chloride-Dieihyl ether at 298.0 K temperature. 


;h) 



M, ETHANOL — - ETHANOL 
PRESSURE OF THE SYSTEM (bar) 


1.0130 '■ 

CHANGED TSONOPOULOS-K VALUES 


.0000 

.0150 







.0150 

.0000 







tempCc) 

X 


Y. 

exp 

Y. 

cal 

t (exp) 

t (cal) 

xl 

x2 

yi 

y2 

yi 

y2 

77.0000 

77.0000 

.0730 

.9270 

.1200 

.8800 

. 1121 

.8879 

75.8000 

75.8000 

.1420 

.8580 

.2220 

.7780 

.2108 

.7892 

74.0000 

74.0000 

.2480 

.7520 

.3620 

-'.6380 

T3500 

.6500 

72.3000 

72.3000 

.3750 

.6250 

.5050 

.4950 

.4987 

.5013 

70.9000 

70.9000 

.4700 

.5300 

.6000 

.4000 

.5986 

.4014 

69.2000 

69.2000 

.6000 

.4000 

.7250 

.2750 

.7204 

.2796 

67.6000 

67.6000 

.7250 

.2750 

.8200 

. 1800 

.8226 

. 1774 

65.8000 

66.0000 

.8880 

.1120 

.9320 

.0680 

.9352 

.0648 


AVERAGE DEVTATTON TM TEMPERATURE Cc) 
.025 

AVERAGE DEVIATION IN Y-VALUES, 

DYl DY2 

.OOGIG .OOGIG 


T'ahlc 5.1 


15: ()()iiipa,iiH()ii ol’ cxpcriiucutal ami calculated values ol ]){ ami p 


Wong-Sandlcr mixing rules for Methanol- Ethanol at 1.013 bar pressure 


• 11 ) 


from 



ACETONE — CHLOROFORM — HEXANE 
TEMPERATURE OF THE SYSTEM (K) 


308.0000 

CHANGED TSGNQPOUI.OS-K VALUES 


.0000 

0000 .1500 






.0000 

0000 .1000 






.1500 

,1000 .0000 






PRES.. 

. (mm) 


X 

Y. 

. exp 

Y. 

..cal 

P(exp) 

P(cal) 

xl 

, x2 

yi ■ 

y2 

yi 

y2 

300.5000 

319.4629 

.0650 

.6000 

.0600 

.6000 

.0643 

.6022 

320.9000 

315.3640 

.0870 

.1000 

.2670 

.0890 

.2667 

.1277 

287.0000 

305.3841 

.1090 

.7990 

.0600 

.7930 

.0779 

.7771 

280.1000 

301.0839 

.2070 

.7000 

.1580 

.6790 

.1768 

.6570 

275.6000 

306.1683 

.3050 

.6010 

.2760 

.5450 

.2894 

.5222 

344.3000 

361.7745 

.4000 

.3040 

.4620 

.1840 

.4512 

. 1985 

340.7000 

360.8897 

.4980 

.3050 

.5410 

.1730 

.5119 

.1878 

413.9000 

411 . 1522 

.5990 

. 1030 

.6210 

.0420 

.6081 

.0491 

401.7000 

404.5984 

.7100 

. 1070 

.6690 

.0430 

.6553 

.0500 

383.2000 

383.8397' 

.8070 

. 1040 

.7520 

.0440 

.7401 

.0516 

AVERAGE DEVIATION IN 

PRESSURE (mm Hg) 






13.83855 

AVERAGE DEVIATION IN Y-VALUES 

DYl DY2 DY3 

.01331 .01575 .00940 

Table 5 1 Kr Coinparisoii of experimental and calculated values of y.- and p from 
Wong-Sandlcr mixing rules for Acetone-Chloroform-Hexane at 308.0 K temperature. 
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ACF/rnNE --- OIII.nRGFnRM — - TOLUENE 
PRESSURE OF THE SYSTEM (bar) 


1.0130 

CHANGED TSONOPOULOS-K VALUES 


.0000 

.1000 

.1700 






.1000 

.0000 

.0000 






.1700 

.0000 

.0000 






TEMP 

• CO 

X 


Y. 

. exp 

7. 

cal 

t (exp) 

t (cal) 

xl 

x2 

yi 

y2 

yi 

y2 

88.8000 

84.5000 

. 1.000 

.1000 

.4130 

. 1720 

.4030 

.2224 

78.4000 

74.0000 

.2000 

.2000 

.5180 

.2670 

.4980 

.2957 

69.4000 

64.5000 

.2000 

.6000 

.2650 

.7000 

.2856 

.6546 

70.3000 

65.0000 

.3000 ' 

.4000 

.4860 

.4330 

.4725 

.4370 

OR.nooo 

6-1 . 0000 

,^000 

. 3000 

.6270 

.2850 

.6019 

- . SOTO- 

66.6000 

62.5000 

.5000 

.2000 

.7660 

. 1600 

.7149 

. 1943 

64.8000 

62.0000 

.6000 

. 1000 

.8540 

.0640 

.8144 

.0929 

62.4000 

60.0000 

.7000 

.1000 

.8880 

.0610 

.8430 

.0897 

61.8000 

59.5000 

.7000 

.2000 

.8450 

. 1320 

■“.“7849 

.1793 

60.4000 

58.5000 

.8000 

.1000 

.9170 

.0600 

.8732 

.0889 


AVERAGE DEVIATION IN TEMPERATURE CO 
3.690 

AVERAGE DEVIATION IN Y-VALUES 

DYl DY2 DY3 

.03291 .03195 .01571 


I’ahUi 5.1.17: (lompariHoii of cxpcriinciiUl and 
Wonfi-Sandlr.r mixing rules for Acctone-CMlorojo 


calculaicd values of yi and t from 
rm-Tolucnc ai 1.013 bar pressure.' 



ACETUNl:; --- 'i'E'i'HACHLUHDMETHANE BENZENE 
PRESSURE OF THE SYSTEM' (bar) 

1.0130 

CHANGED TSQNGPUIJLOS-K VALUES 


.0000 

.0000 

.0500 






.0000 

.0000 

.0500 






.0500 

.0500 

.0000 






TEMP 

• CO 

X 


Y. 

exp 

Y, 

. cal 

t (exp) 

t (cal) 

xl 

x2 

yi 


— yi 

y2 

79.7000 

79.5000 

.0160 

.0170 

.0523 

.0200 

.0516 

.0253 

75.3000 

74.2500 

.0520 

.2630 

.1450 

.2580 

.1376 

.2917 

65.8000 

68.0000 

.1410 

.7990 

.3460 

.6140 

.3204 

.6275 

65.1000 

65.0000 

.2880 

.2860 

.4800 

.2100 

.4748 

.2415 

64.2000 

64.5000 

.3130 

.2440 

.5240 

.2100 

.5009 

. 2082 

63.0000 

62.5000 

.4700 

.0910 

.6320 

.0830 

.6424 

.0801 

60.1000 

60.2000 

.5600 

.1960 

.6740 

. 1760 

.6856 

. 1584 

58 . 5000 

59 . 5000 

.B-lOO 

.4310 

.6820 

.3030 

.6653 

.3151 

58 . .3000 

59.0000 

.6310 

.2960 

.7320 

.2280 

.7243 

.2273 

57.5000 

58.0000 

.7420 

.1770 

.7950 

. 1600 

.8017 

. 1440 

57.4000 

57.4000 

.8570 

.0630 

.8860 

.0620 

.8899 

.0550 

57.4000 

57.4000 

.7840 

.1860 

.8260 

.1260 

.8280 

.1515 


AVERAGE DEVIATION IN TEMPERATURE (“C) 


.554 

AVERAGE DEVIATION IN Y-VALUES 

DYl S DY2 DY3 

.00839 .01378 .01155 

'I'alilc 5 . 1 . 18 : ( loinparisoii of <:x[H:rinu:iiUl and calculatod values of y,- and t from 
S (iV'd l('V inixin^ rules loi' I civcK'Jilovoitict I kihc- I dcjizc'tic at l.Uld bar 

})ressure. 



MKTIlANOh UKrTANli TULUENE 
PRESSURE OF THE SYSTEM (bar) 


1.0130 

CHANGED TSONOPUULOS-K VALUES 


.0000 

.3000 

.1500 






.3000 

.0000 

.0000 






.1500 

.0000 

.0000 






TEMP. CC) 

X 

.y 

Y. 

exp 

Y. 

cal 

t (exp) 

t (cal) 

xl 

x2 

yi 

y2 

yi 

y2 

62.9600 

64.0000 

.3172 

.2533 

.7458 

.1826 

.7665 

.1199 

63.3400 

64.5000 

.3711 

.0962 

.7797 

.0635 

.7888 

.0615 

59.9700 

59.5000 

.5409 

.3412 

.7520 

.2043 

.7973 

.1692 

61.3500 

60.5000 

.6038 

.1717 

.7697 

. 1412 

.7926 

. 1318 

62.5100 

62.5100 

.6114 

.0911 

.7864 

.0911 

.8048 

.0867 

61.5900 

61.5900 

.7494 

.0787 

.7927 

.1123 

.7956 

. 1183 

59.9600 

59.5000 

.7733 

.1510 

.7563 

.1969 

.7647 

.1986 

60.7700 

•' 60.0000 

.8486 

.0684 

.7815 

.1517 

.7824 

.1604 

AVERAGE 

DEVIATION 

IN TEMPERATURE CC) 






0.596 

AVERAGE DRVIAT.TOM IN Y-VALUES 

DYl DY2 DY3 

.01609 .016252 ‘ .01443 

Table 5.1.19: (Jom[)aiison of cxperiinciilal and calculalcd values of yi and t fipra 
Wong-Sandkr mixing rules for McLiianol- Hep Lane-Toluene at 1.01 3 bar pressure. 


-I'l 



MKI'IIANDl, rKTHAOllUJlUJMia'llANl!: UENZliNK 

TEMPERATURE OF THE SYSTEM (K) 

307.6800 

CHANGED TSONOPOULOS-K VALUES 


.0000 

.1500 .1500 







.1500 

.0000 .0000 







.1500 

.0000 .0000 







PRESS 

. (mm) 

X 


Y. 

,exp 

Y. 

. cal 

P(exp) 

P(cal) 

xl 

x2 

yi 

y2 

yi 

y2 

291.1100 

284.8949 

.2075 

. 1900 

.4920 

. 1472 

.5077 

. 1538 

302. 1300 

294.5722 

.2110 

.3879 

.4804 

.2774 

.4829 

.2911 

308.6300 

300.1083 ' , 

. 1987 

.5876 

.4733 

.3999 

.4681 

.4058 

307.2300 

307.4335 

,3781 

.3122 

.5043 

.2747 

.5284 

.2697 

308.1300 

308.2517 

,5543 

.2078 

.5308 

.2501 ' 

.5512 

.2425 

298.8000 

302.2973 

,7599 

. 1076 

.5903 

.2080 

.5945 

.2145 


AVERAGE DEVIATION IN PRESSURE(mm Hg) 
4.35286 


AVERAGE DEVIATION IN Y-VALUES 

DYI DY2 DY3 ^ 

.01201 .00755 .01115 

'I'ahlc^ .5.1.20: (lompiulKoii of (;x[)erimonl.al and calculated values oi y,- and p from 
Wo'iifj' mixing rnl('H loi' Mvlhaiiol-'rclvnclilovovK'lhtinv-lh'iizvnv. at UOi.dS K 

tempera, ture. 
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Temperature (^C) — Temperature (oC) 


* 


vs. ?/ 


Mole fraction, iji — ^ 

Figure 5.2.1; Comparison of experimental and calculated values of yi and t from 
Wong-Sandler mixing rules for Acetone-Benzene at 1.013 bar pressure. 


vs. 


vs. 


0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

Mole fraction, i/i — *• 


Figure 5.2.2: Comparison of experimental and calculated values of y,- and t from 
Wong-Sandler mixing rules for Acetone-Chloroform at 1.013 bar pressure. 








Mole fraction, 'iji 


Figure 5.2.5; Comparison of experimental and calculated values of yi and t from 
Won(j-S(indler mixing rules for Acctone-Tetrachloromcthane at 1.013 bar pressure. 



Figure 5.2.6: Comparison of experimental and calculated values of ?/,• and t from 
fFonY/-.SVind/f:r mixing rule.s for darbon disulfide- Acetone at 1.013 bar pressure. 









0 0.2 0.4 0.6 0.8 1 


Mole fraction, iji — > 

l*'igure 5.2.9: Comparison of experimental and calculated values of yi and t from 
Wo'tKj-S (i.'ndlcv ix\\y.\ng rules lor A cct-aldchydc- Benzene i\l l.Old ba.r pressure. 



Figure 5.2.10: Comparison of experimental and calculated values of ?/,■ and t from 
Worifj-Sandler mixing rules for Diethyl elher-Benzene at 1.016 bar pressure. 






l'’igurc 5.2.11: Conipaxison of experimental and calculated values of yi and p from 
Wong- Sandler mixing rides for Diethyl ctlicr-ChloroJorin at 290.0 K temperature. 



Figure 5.2.12: Comparison of experimental and calculated values of y,- and p from 
Wong-Sandler mixing rules for Diethyl ether-Carbon disulfide at 293.0 K tempera- 

ture. 


CENTRAL LtBRAR> 

KAMniR 





rressure (ram n 


Mole fraction, y\ 


I'Mgnre 5.2.13: Comparison of experimental and calculated values of y,- and p from 
Wong-Sandlcr mixing rulc.s for Diethyl cihcr-Elhyl ace.Lalc at 303.0 K temperature. 



Pdgure 5.2.14: Comparison of experimental and calculated values of yi and p from 
Wong-Sandler mixing rules for Ethyl chloride-Diethyl ether at 298.0 K temperature. 
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I'Mgure 5.2.15: Comparison of ex 
Woncj- Sandler mixing miles for M 
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I'^igure 5.2.16: Comparison of ex 
Wona-Sandler mixine: rules for Ai 





lemperature 



Figure 5.2.17; (!oiu[)ari.soii of experimental and calculated values of m and t from 
Woruj-Sandla- mixing rules for Acelone-Cliloroform-Tolucne at 1.013 bar pressure. 



Figure 5.2.18: Comparison of experimental and calculated values of yi and t from 
Wong-Sandler mixing rules for Acetone- TcLrachloromethane-Benzerie at 1.013 bar 

])ressure. . , . . . 





rressiire [mm n 



' if^urc' 5.2.19; ('onipiirisoii 
Wong-Sandlcr mixing rules 


of experimental and calculated values of y.- and t from 
for A'lelhanol-lIeptane-Toluene at 1.013 bar pressure. 



rure 5.2.20: Comparison of experimental and calculated values of yi and p from 
Ing-Sandler mixing rules for Mclhaml-Tetrachloromethane-Bmztne 307.68 

noerature. 





Chapter 6 

CONCLUSIONS 


The equations of state which are applicable to both vapour and liquid phases 
have obvious advantages over activity coefficient models especially when one is inter- 
ested in large ranges of temperature and pressure including critical and supercritical 
regions. However phase equilibrium calculation using equations of state are rather 
lengthy except for simple equations of state. Another problem with the use of the 
equations of state technique is the insufficiency of information about the exact form 
of equations of state for mixtures. 

While calculating the cross second virial coefficient by Tsonopoulos correlation, 
binary interaction parameters for cross temperatures have been tuned in order to 
match calculated results with experimental results. This is one of the disadvantages 
of this model. But one can overcome this difficulty by using optimization technique 
(e.^., genetic algorithm, direct search methods, gradient-based methods) to determine 
binary interaction parameters for cross temperatures and these values can be used 
for other systems at different conditions. 


56 



Bibliography 


1. Eubank, P. T., Shyu, G. and Hanif, N. S. M., New procedures for appli- 
cation oj the Wong-Sadler mixing rules to the prediction of vapour- liquid 
equilibria, Ind. Eng. Chem. Res., 34, 314-323, (1995). 

2. Gmehling, J., Onlien, U. and Arlt, W., ’’Vapour Liquid Equilibrium Data 
Collection” Vol. 1, parts 2a, 2b, 34-4. 

3. Huron, M. J. and Vidal, J., New mixing rules in simple equations of 
state for representing vapour-liquid equilibria of strongly nonideal mixtures, 
Fluid Phase equili., 3, 255, (1977). 

4. Orbcy, II. and Sandler, S. I., On the combination of equation of state and 
excess free energy models, Fluid Phase equili., Ill, 53-70, (1995). 

5. Panagiotopoulos, A. Z. and Reid, R. C., New mixing rules for cubic equatins 
of state for highly polar asymmetric mixtures, ACS Symp. Ser., 300: 571- 
582, (1986). 

6. Peng, D. Y. and Robinson, D. B., A new two-constant equation of state, 
Ind. Eng. Chem. Fundam. 15, 59, (1976). 

7. Prausuitz, J. M., Lichtenthaler, R. N. and de Azevedo, E. G., ’’Molecular 
Thermodynamics of Fluid Phase Equilibria”, 2nd Ed., Prentice-IIall Inc., 
1986. 

8. Reid, R. C., Prausuitz, J. M. and Poling, B. E., ’’The Properties of Gses 
and Liquids”, 4th Ed., McGraw-Hill International Editions, 1987. 

9. Schwartzentruber, J. and Renon, H., Extension of UNIFAC to high pres- 
sures and temperatures by the use of a cubic equation of state, Ind. Eng. 
Chem. Res., 28: 1049-1055, (1989). 

10. Smith, J. M. and Van Ness, H. C., ’’Intruduction to Chemical Engineering 
Thermodynamics”, 4th Ed., McGraw-Hill, New- York, 1987. 

11. Soave, G., Equilibrium constants from a modified Redlich-Kwong equation 
of state, Chem. Eng. Sci., 27, 1197, (1972). 

12. Stryjeh, R. and Vera, J. H., PRSV: An Improved Peng-Rohinsosn equation 
of state for pure compounds and mixtures. Can. Jl. Chem. Eng., 64: 
323-333, (1986). 

13. Tsonopoulos, C., An emperical correlation of secod virial coefficients, AIChE 
Journal, 20, 2, 263, (1974). 


57 



14. Wong, D. S. H., Orbey, H. and Sandler, S. 1., Equation of state mixing rule 
for nonideal mixtures xising available activity coefficient m,odel parameters 
and that allows extrapolation over large ranges of tempe.r'aturc a.nd pressure^ 
Ind. Eng. Chern. Res., 31, 2033-2039, (1992). 

15. Wong D. S. H. and Sandler, S. I., A thermodynamically correct mixing rule 
for cubic equations of state, AlChE Journal, 38, 5, (1992). 


58 



